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CRITICAL LEVITATION LOCH FOR SPHERES ON CRYOGENIC FLUIDS 


by R. C„ Hendricks and S. A, Ohm 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

The forces and governing equations for spheres floating on a fluid of 

lower specific gravity are examined. Three similarity parameters, Bond 

9 

number (Bo - (p - p„)gR“/ g_), component density ratio of the sphere to 

i v O C 

the fluid (p_ - p/p, ~ Ptt)j wetting angle (a) arid the critical levitation 
locii, established herein, are sufficient to determine whether or not the 
sphere will float. Data appear in good agreement with the analysis for a 
limited range of Bond numbers and a wetting angle of it. These data and 
the theoretical curve for a ~ n may be approximated by 

^IL-Ox = 1,05 +i^i§ 0. 1 < Bo < 30 

P t - P v Bo 

Those spheres with density ratios less than this criteria will float , others 
will sink. More data are needed to verify the .analysis over a larger range 
of Bond numbers and wetting angles, 

INTRODUCTION 

Floating a waxed steel pin on water has long been a source of amuse- 
ment and an apparent paradox to Archemedes principle. However, a 
closer examination reveals that one simply has not accounted for all the 
forces acting to float the pin, A large contribution comes from the effects 
of surface tension which acts in two ways: first, the free liquid- vapor 
interface is curved above the solid-liquid-vapor contact line; as such it 
acts much like a retaining wall, and holds back a fluid head; second, the 
solid-liquid interfacial, tension forces acting on the body are transmitted 
to the solid- liquid- vapor contact line and held in equilibrium by the free 
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liquid- vapor interface . An analogy for the latter forces is that of a thin 
shell tank being supported at the bolting flange. 

Beside the usual curiosity of the phenomenon, floating of a noncry- 
ogenic fluid on a cryogenic fluid (and vice-versa) is of significant impor- 
tance in fuel spills (ref. 1). Floatation phenomena are also of interest in 
studies on controlling hypergolic fires (ref. 2); in the preservation (or 
destruction) of biological specimens, such as in blood preservation and 
in enzyme preparation (ref. 3). Floatation phenomena could be of signif- 
icance in the production of monolayer films, pollution control, and as 
fundamental tools for the determination of surface wetting characteristics. 

The shapes of axisymmetric interfaces have been calculated by Huh 
and Scriven (ref. 4) and Hendricks and Baumeister (ref. 5)„ 

Huh and Scriven (ref. 4) numerically solved the Laplace capilliary 
equation, analyzed the errors involved, and presented solution charts and 
tables for the interface configurations. Hendricks and Baumeister (ref. 5) 
determined the heat transfer characteristics of a water-sphere floating on 
a sea of liquid nitrogen in Leidenfrost film boiling. Nutt (ref. 6) examined 
the forces required to float an isothermal sphere. However, to approxi- 
mate the submergence depth, Nutt (ref. 6) neglected the curvature along 
the line of contact and used the closed form solution for a cylinder. Huh 
and Scriven (ref. 4) found curvature along the line of contact to be impor- 
tant for axisymmetric surfaces at low values of Bond number. While ref- 
erences 4, 5, and 6 deal with the problem of floating spheres and calcu- 
lating interfaces, only reference 6 gives an approximate criterion for 
determining if a given isothermal sphere will float on a specified fluid. 
Such a criterion has not been established for nonisothermal floating 
spheres and a more accurate solution needs to be established for the iso- 
thermal floating spheres. In the ensuing sections we will formulate the 
governing equations, determine the similarity parameters involved, 
numerically determine the optimum levitation locii (i.e. , a float or no- 
float criteria) and examine the asymptotic behavior of the solution for 
both the nonisothermal and isothermal floating spheres. Some experi- 
mental data are presented to check the analysis. 
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THEORETICAL ANALYSIS 

In this section we wish to determine the conditions permitting a fluid 
to support a sphere having a higher specific gravity than its own. The fol- 
lowing analysis will produce three basic parameters which together define 
the maximum floating conditions: Bond number (Bo), that is, the ratio of 
buoyancy to surface tension forces; the wetting angle of the interface (a)*; 
and the ratio of solid-liquid specific gravities (p - p /p - p Tr ). 

We proceed with analysis of the forces at the sphere-liquid interface 
in order to determine the optimum levitation loci; that is, the conditions 
permitting floatation of a maximum density sphere. We consider first the 
force balance on the nonisothermal floating sphere in film boiling (fig. 1(a)), 
then the force balance on an isothermal sphere (fig, 1(b)), incorporating the 
work of references 4 and 6. Both cases are subject to the constraint of min- 
imum energy conditions at the interface, A consideration of the limiting 
cases then defines the asymptotic solutions of the desired loci, 

Basic Force Balance for Nonisothermal Sphere 

As shown in figure 1(a), the basic model for the nonisothermal floating 
sphere in film boiling includes a vapor envelope separating the floating sphere 
and the supporting liquid. This situation occurs, for example, when a water 
sphere at room temperature (295 K) is placed on saturated liquid nitrogen 
(77.4 K), and floats on a cushioning layer of nitrogen vapor. 

The balance of forces acting to support the sphere, see figure 1(a), may 
be expressed as: 


tThe value of wetting angle for an advancing interface and that of a re- 
ceding interface is herein referred to as a, i,e. , no distinction is made in 
the analysis. Thus when using the results, one must use the proper wetting 
angle; namely, advancing or receding. Also, no allowance for the variation 
in wetting angle with time is made, although experiments with time dependent 
values of a (or surface tension a, or both) could still use the results pre- 
sented herein by considering time increments. 
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The pressure in the vapor gap at the sphere surface, P , may be expressed 
in terms of the pressure drop across the vapor gap, APg , the pressure 
drop across the interface, AP^, and the liquid pressure acting on the inter- 
face P^ , see figure 2(a) 


P — P. + AP - APt. 

Ol CTO 


Using figure 2, the expressions for P and P become 

6 ” 
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Hf + Z o - R q (cos co + cos 0*) 


H lM + [ Z o 


(R rt + 6) (cos 9 - cos 
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Pyg 


0*)] 
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( 3 ) 

(4) 

(5) 


and because the interface is nearly spherical, the pressure drop across the 
interface becomes 


AP = 2ct (6) 

a (R 0 + 6) 

The pressure drop across the vapor gap AP § is determined from equa- 
tions of reference 5, and equation (4), 


fThe difference in vapor head is quite small and was omitted in ref. 5. 
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AP 5 = R o^ " p 



R o (P; - Py)g W s 

ffW*f(0*) 


X 




+ Z q - R q cos 0*) - R cos 
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( 7 ) 


where 


<P = 


1 + 6 



(8) 


The shear t ^ can be evaluated using the velocity profiles and non- 
dimensional form of reference 5„ 

J r r0 dA 3 = (1 - q>) W s (9) 

Substituting equations (5), (6), and (7) into equation (3) and then sub- 
stituting equations (2), (3), (4), and (9) into equation (1) gives 


0 = q)W c 


+ 2 n 


■ 0 * 


9 

(H^ + Z Q - R q cos 0*) - R q cos cc|p v gR^ cos u> sin c o do> 
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o 
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w)/ 



(R + 6) cos 0* 
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P,g 


R p(P; - Py)g w s /l - 
irw*f(0*) V cp 


■(Hj + Z Q - R o cos 0*)p y g(cp - 1) 


+ PoP v (f " 1) + <R 0 + 6)P^ 


g cos 0) R^ cos 9 sin 9 d0 


( 10 ) 


If the parameters Z Q and 9* were known, then equation (10) could be 
integrated in a straight forward manner, (see appendix B) and the force 
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balance would be complete. However, these parameters must be deter- 
mined by a simultaneous solution of this force balance and the interface 
equations, to be found in the theoretical section - Structure of the Interface. 

Basic Force Balance for the Isothermal Sphere 

As is shown in figure 1(b), the forces in the isothermal sphere-liquid 
system act di rectly on the floating sphere. The force balance, equating 
the weight of the sphere with the buoyancy and surface tension forces, is 
given by a static force balance analogous to equation (1): 

W s - F B - F o = 0 <“> 

Again, the weight of the sphere is, 


w < 12 > 

However, the surface tension forces acting on the sphere are considered to 
act along the contact circle giving use to the surface tension force, see 
appendix B, 

F a = - f cr^ y sin (0* + a) dl 

- -2ttR 0 ct^ v sin 0* sin (0* + a) (13) 

where dl = R q sin 0* d<& over the contact circle. „ A similar argument is 
given by Nutt, reference 6. 

Appendix C includes a discussion and definition of this surface tension 
force in terms of adhesive surface tension. The buoyancy force of equa- 
tion (11) is given by the difference between the supporting liquid force on 
the liquid-sphere area A*, and the downward vapor pressure on the vapor- 

sphere area A^. 



where and P^ are defined by equations (4) and (5), see figure 2(b). 
Again, the buoyancy contribution of equation (14) to the force balance of 
equation (11) is a matter of simple integration, given the parameters Z Q 
and 0*. We now investigate the interface equations in order to solve for 
these parameters and simplify the force balance in the theoretical section - 
Force Balance for Both Isothermal and Nonisothermal Levitated Spheres. 

Structure of the Interface 

In determining the minimum energy configuration of the liquid-vapor 

interface in both the isothermal and nonisothermal cases, we investigate 

the forces at the interface, which is considered to be a thin shell as seen 

in the surface element of figure 3, enlarged in detail AA.t For values of 

0 greater than 0*, the interface seeks a minimum energy configuration 

subject to two constraints: (a) the balance of forces on the sphere at the 

load circle, and (b) the boundary conditions at °o. Referring to figure 3, 

2 

the balance of forces in the X -direction becomes: 

3r«n_ 

— — — dT d| - (R.Nj- d£ dT) cos (ir - r) + dr d| = 0 (15) 

dr 1 * * 

3 

and in the X -direction: 

-R 3 N t d£ dr + (R-jN| dT d£) sin (n - t) + (P^ - P v )RgR^ d| dr = 0 (16) 

where N r and are surface forces per unit length in the £ and t 
directions, respectively. It can be shown that these equations, (15) and 
(16), represent the governing equations of the liquid interface if 


"fin ref. 4, Huh and Scriven consider various forms of the interface equa- 
tions and provide an excellent review of previous works . 
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n s = n t = v ( 17a ) 

and 

Y f = 0 (17b) 

See references (5) and (8). These governing equations may be written as 
two first-order, nonlinear, ordinary differential equations, subject to 
equation (17), and reduced to; t 

D x u = -- - z (18) 


where 


D x z = 


u 




= tan t 


P ( - P v = Z ( f>l - p v )g 


z 

Z = — = 


a 


Z v 


(P; - p v )e 


(19) 


X = 


X 


u = sin r 


Now, equations (18) and (19) are subject to the boundary conditions ; 
Sphere load circle. - 


x = x Q = i\Jbv) sin 0 ! 


z = z. 


o 


( 20 ) 

( 21 ) 


'See refs. 4, 5, and 8. 
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t o = (n - 9*) + (tt - a) (22) 

Interface at <*> . - As x — °° (L e. , x = M » 0) 

z 0 (23) 

r — it (24) 

where 

R n 

Bo = — = Bond number (25) 

L 2 

The (-) sign in equation (19) is used when r < tt/2 and (+) sign used 
when r > ij/ 2. As the conditions at t = u/2 represents a singularity, this 
region was handled by changing the independent variable from x to t„ 
When 0.93 < u ^ 1, equations (18) and (19) were rewritten as 


D t x 


D t z = - 


X COS T 

xz + sin r 

(18a) 

x sin t 

(19a) 

xz + sin r 



Details of a more sophisticated and accurate technique are found in refer- 
ence 4. In either case, a fourth-order Runga-Kutta technique forms the 
basis of the numerical solution. 

Force Balance for Both Isothermal and Nonisothermal Levitated Spheres 

For a given Bo and a, the interface equations (18), (19), (18a), and 
(19a) yield solutions for z q and 0*, which can then be used to evaluate 
either equation (10) or equation (14), thus completing the solutions of equa- 
tions (1) and (11). 

Equation (10) is evaluated in appendix B in a more precise manner 
than reference 5. For the nonisothermal floater, 5/R Q « 1 which can be 
demonstrated using the results of reference 5. This implies that 
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and equation (BIO) reduces to an expression in terms of p, Bo, and 9* 





9 * 


- — cos 0* sin^ 
2 


9 * 


(26) 


Integration of equation (14), gives for the buoyancy force on the iso- 
thermally floating sphere: 




pj + j\ 

3( f>l - P T ) 


?° B J ± + £° ? 3 J * + ?o_(i . cos 2 0*) 
2 6 2R 0 


(27) 


Substituting equation (13) for surface tension force and equation (27) 
for buoyancy force into equation (11) gives the force balance on the iso- 
thermal sphere, expressed in terms of p, Bo, 9 *, and 

2 fe . '. fy = 1- I COS 0* +1 cos 3 0* + Z ° Sln - 2 .— - sln 9 * sln < a + 9 *> 

3 p^-p v 3 2 6 2 a/Bo Bo 

(28) 

Note that equation (26) for the nonisothermal floating sphere is iden- 
tical to equation (28) for isothermal floater, for a - ir, assuming of 
course that 8 /R q « 1. Thus for the nonisothermal case (film boiling), 
we can validly employ the maximum floating criteria derived for the iso- 
thermal case, with a = n. 


tThe form of eq. (28) can be altered to agree with that of Huh and Scriven, 
ref. 4, and is similar to the results of Nutt, ref. 6, within the framework of 
his assumptions. 
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Critical Levitation Loci and Asymptotic Limits 

In order to obtain graphs of levitation loci as a function of density ratio 

(p_ + p /p, - p._) over a range of Bond numbers (Bo) for given values of wet- 

ting angle ( a ), we solve equations (18), (19), (18a), and (19a) subject to the 

boundary conditions, equations (20) through (24) and the constraint equation 

(28) . The solution is obtained by specifying the Bond number, wetting angle, 

and the 9* and z parameters,, Once a fourth order Runga-Kutta solution 

u 

is attained, the above procedure is repeated for a range of values of 9* and 
z Q until the ratio (p g + p y /p^ - p y ) of equation (28) reaches maximum value. 
Figure 4(a) presents resulting maximum density ratio loci over Bond number, 
Bo, values for a = it, tt/ 2, and 3n/5 (3rr/5 is for Teflon floating on water); 
figure 4(b) gives corresponding submergence angles 9*, and figure 4(c) the 
corresponding submergence depths z , as functions of Bo. For a given 
Bo and a, spheres with density ratios less than specified by this locus will 
float; however, spheres with density ratios greater than specified by the 
locus will not float. 

Investigation of the characteristic trends of an equation is both instruc- 
tive and expedient in determining the desired results. We now investigate 
the two limiting cases of very large and very small Bond number limits. 

As Bond number becomes very large (Bo — °°) equation (28) is approxi- 
mated by : 


' P s "-P V ' 
3 Vp* - Py/ 


i - i cos 8* + A cos 3 8* 
3 2 6 


As Bond number becomes quite small, (Bo — e « 0), (i.e. , 
forces dominate), equation (28) becomes 

2 f p s f Py\ _ z o sin ® ' sin 9 * sin (9* + a) 
2 /\/Bo Bo 

sin 9* sin (9* + a) 



(29) 

surface tension 

(30) 


Bo 


(31) 
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The expected asymptotic behavior of the critical Levitation Locii as 
determined from equations (29) and (31) is summarized in table I. 

The asymptotic limits established in this section. are in good agree- 
ment with the computed locii, and establish the trends for the limits of low 
and high Bond numbers (see fig,. 4) , 

RESULTS AND DISCUSSION 

A number of solid and liquid spheres have been floated on water and 
liquid nitrogen, respectively, to test the Optimum Levitation Locii of 
figure 4, The solid teflon spheres and fluid containers were carefully 
cleaned prior to attempting to float a given sphere in the isothermal state. 
The spheres were then placed on the supporting fluid with tweezers. 

While this technique is satisfactory for these experiments, more precisely 
controlled experiments will require a sting or comparable type support in 
conjunction with an analytical balance. In the experiments of Nutt (ref. 6) 
a centrifuge was employed to vary the local acceleration and give first 
order results; however while the technique is a good one, a larger con- 
tainer should be employed to minimize the effects of curvature at the con- 
tainer wall. The liquid spheres were formed using a hypodermic needle, 
squeezed out and carefully dripped onto the liquid nitrogen surface (the 
nonisothermal state). While this technique is adequate for our case, more 
accurate data will require a technique which can minimize the surface 
impact. The size of the frozen spheres were measured by comparison to 
teflon spheres and checked using a scale. 

To determine if a given material would float on a specific fluid, vari- 
ous diameter spheres were tried. Sphere diameter was increased until 
the surface could no longer support the sphere and it sank. In most cases, 
the graduations in sphere diameter were too coarse to permit the determin- 
ation of the critical size. These "largest size" levitated sphere data are 
presented in table U and figure 4(a), 

The data for a - n (Le. , those spheres floating on liquid nitrogen) 
are in good agreement with the theoretical line for a limited range of Bond 
number. The teflon sphere data, a ~ 3n/5, also appear to be in good 
agreement with the analysis. 
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To illustrate the sensitivity of the data, m 0. 31 centimeter diameter 
carbon tetrachloride sphere would float momentarily on liquid nitrogen and 
sink to the bottom while a 0,28 centimeter sphere floated. Carbon tetra- 
chloride spheres 0.31 centimeter diameter and above are classified as non- 
floaters on liquid nitrogen. A 0.635 centimeter diameter teflon sphere 
would not float on water while other teflon spheres (0. 318, 0. 3S7, and 0,478 
cm diam. ) would float on water. 

From these results it appears that the parameters of Bond number 
(Bo = (jp z - p v )gR|j/ag c ) ? density ratio (p g - p y /p z - p y ), and the wetting 
angle (a) used in conjunction with figure 4 are sufficient to predict whether 
or not a given sphere will float on a specified fluid. A decrease in the wet- 
ting angle (increased wetting which effects the surface tension support) 
while holding Bond number constant will affect a decrease in the density 
ratio (p g - P v /p z - P v ); this means that a less dense sphere can be sup- 
ported. It is apparent from figure 4 and the data that the effects of a 
become more significant at low Bond numbers, and a contaminated surface 
will greatly alter the results. As a more specific example, teflon spheres 
of diameter greater than 0. 318 centimeter would not float as received from 
the manufacturer. However, degreasing and cleaning the surface resulted 
in good agreement with the optimum levitation locus for teflon on water. 

The effects of advancing and receding wetting angle, a, have been in- 
corporated in the analysis because the a-parametric curves of figures 4(a) 
and (b) represent either value of a. The time dependency of a has not 
been incorporated; however, If the values of a are known as a function 
of time, figures 4(a) and (b) should give good results. The solutions should 
also be valid for temperature variations in surface tension cx Zy of the sup- 
porting fluid provided the heating rates are not large or the fluid near its 
thermodynamic critical point (Le. , the sphere and the liquid remain close 
to isothermal and heat transfer is very small). If this premise is violated 
then the. governing equations must include heat transfer. 

Bond number appears to be the most easily controlled parameter and, 
as it characterizes the supporting liquid interface, becomes the fundamental 
independent variable in the experiments. (e. g. , the size of the sphere was 
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changed in our experiments and acceleration was the variable in the experi- 
ments of Nutt (ref. 6)); the value of o Jv could be changed by heating the 
fluid. Where the Bond number is changed by changing sphere size, the 
values of 0* and z Q can be estimated by the technique presented in appen- 
dix D. Other fluids and solids need to be examined to determine the range 
of validity of the analysis. 

CONCLUSIONS 

The forces and associated governing equations for isothermal and non- 
isothermal spheres (film boiling) floating on a fluid of lower specific gravity 
have been examined. Criteris for floating spheres called the Critical Levi- 
tation Locii have been established. Bond number Bo = \jfl ” p v^ gR ol/ /ug c ; 
density ratio, (p_ - p /p 7 - pj; wetting angle, a; used in conjunction with 
figure 4 appear to be sufficient to determine if a given sphere will float on 
a specified fluid. For the case of a = tt the Critical Levitation Locus may 
be approximated by: 

0l = 1.05+ii£5 0.1 < Bo < 30 

P t - P v Bo 

Data for glycerine, carbon tetrachloride, and water spheres floating 
on a sea of liquid nitrogen appear to be in good agreement with the analysis. 
Data, for teflon spheres on water also appear to be in agreement with the 
analytical results. 

More exact and extensive measurements are required to determine the 
validity of the Critical Levitation Locii over a range of Bond numbers (Bo) 
and wetting angles (a). 

The effects of wetting angle, surface tension, and their time, temper- 
ature dependencies bear further investigation. 
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APPENDIX A 
SYMBOLS 

A 2 

A area, cm 

A q parametric form defined by eq. (B2) 

A 1 (r 3 + (aRg/ar) d^ d£ [N t + (9N t / 3t) dr] , see fig. 3 

B parametric form defined by eq. (B3) 

2 

Bo Bond number, (p z - p y )gR 0 /ag c 

Fg lift force due to buoyancy, dyne 

F lift force due to surface tension, dyne 

o - 

f(0*) functional form in eqs. (7) and (10),, cos^ 9 * - 2 cos 9 * + 1 

2 

g local acceleration of gravity, cm/sec 

g c gravitational conversion constant, l[(g)(cm)/(dyne)(see^)] 

arbitrary reference pressure head, cm 

hj arbitrary reference dimensionless pressure head, H^/L 

L reference length, yg^/CP; - P v )g cm 

M a large positive number 

N Surface force per unit length, dyne/cm 

2 

P pressure, dimensional; dyne/cm 

2 

P Q pressure acting on the sphere within the vapor gap, dyne/cm 

P pressure acting on the interface from within the vapor gap, 

9 

dyne / cm 

2 

AP^ pressure drop across the liquid vapor interface, dyne/cm 

2 

APg pressure drop across the vapor gap, dyne/cm 

p pressure, dimensionless 

R sphere radius, cm 

o 
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R 


V 




R 


3 


r 


u 

u* 



w 




X 

Y 1 y 2 ° 


X' 


x 

y j 

z 

Z 


Q! 

6 

e 

C 

e 

9 * 

v 


a 

T 


Gaussian radii of curvature, cm 

radius defined in fig. 3(b) 

radial coordinate, cm 

dependent variable of eq. (18), sin r /.; 

reference velocity, n/R Q cm/sec 

weight of the sphere, dyne 

2 2 

reference weight, pu* R /g c dyne 
horizontal coordinate, dimensional; cm 

coordinate system, fig. 3(b) 

dimensionless horizontal coordinate, X/L 

force tangent to the surface, dyne 

pressure head, dimensional; cm 

dimensionless pressure head, Z/L 

wetting angle, radians 

vapor gap thickness, cm 

a small positive number 

dimensionless radial coordinate, r/R Q 

angular coordinate, rad 

angle to which the sphere is submerged, rad 

2 

kinematic viscosity, cm /sec 

3 

density, g/cm 

surface tension, dynes/cm 

angular coordinate, radians 
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(p dimensionless gap coordinate, 1 + 5/R 0 

<& aximuthal coordinate, rad 

| angular coordinate, radians 

Subscripts: 

adh adhesion 

cr critical 

l liquid 

iv liquid- vapor 

s sphere 

si solid-liquid 

sv solid- vapor 

v vapor 

0 initial 

1,2,3 reference conditions 

r, £ in the direction these angular coordinates 
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APPENDIX B 

SOLUTION OF THE FORCE BALANCE 


In this appendix , equation (10) is integrated and reduced to nondimen- 
sional form. 

Assuming that the parameters Z Q and 9* have been determined else- 
where, the antiderivative of equation (10) becomes: 


0 - 


5 A 

2tj 


+ Vv gR o 



71 - 6 * 


0 


R^> v e 



71 - 9 * 
0 



(Bl) 


where 


A = H, + - R cos 9* 

OlOO 


(B2) 


B 5 H lPy g + [Z 0 - < R o + 6) C0S e *\ p l S 


irw*f{6*) 


f 1 - A - W 

i 1+ -f - 1 " 

\ <P J V 

|V R oj J 


A q (B3) 


Evaluating (Bl) gives 
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W =^ ,2 [1 - cos* 6 *' 


+ r: 


2tt 


o 


)(A 0 P v g - B) - R^p v g /cos “ °* + V 


(or 0 +6)p g- p v g(<p 2 - i)r~ 


R 2 (cos 6* - 1 | (B4) 


0 


2n 



+ cosl£* R 3 g ^i + J_ lp 


R 




o> 


" H l^ - Vv* 1 + 


R, 


o, 


+ r cos e*p g (i +AT +h, p g (l +-5-Y - — — — 

R o/ V V R o +5 


+ V z g 


, , 6 \ _ RJp, - P v )gW s (l - 

R. X +— I cos 0*p, g - -2 — l - 1 5 

' R. / ' ?iw*f(8»),p 


(B5) 


Regrouping the terms of equation (B5) 


0 


w s <P R^g <P 


2v 


K? + V + 22^1! R 2 g[p^ - < 2 - A V ] 


r 2 s in^ 0* r 

^___| Hl p v g(l - „ 2 ) - z 0 g(p v ^ 2 - p t ) +J<L 


COS 0*g 


(R 0 + 6)p ; - R^ p v 


Vi^Ws/l ^ 2 

77W*f(0*) \ <p 


(B6) 


Collecting terms of (B6) and dividing by (p^ - p v ) gR^, 



20 


0 


W s? 

I - Pv>g R o 


- £ 


l p l +p vV\ 

, COS 3 0* 

P Z <P - (2 - <p 2 )p y 

<T 

-.1 

ID 

< 

3 

_ P Z " p v J 


H 


i Pv (ivVm 2 

R o p z " p v \ 2 1 


9* 



l \ dJZ Q* 


sin 


g sin 2 9* + <p^ 

(p t - p v )s «*?> 


cos 0* sin 2 0* / p z ~ Py^ 
2 \ P Z “ p v 


, W s^ - > sin 2 0* 

7TW*f(0*) 2 (p 

Introducing Bond number 


(B7) 


Bo =- 


. K - Py>g R o _ R o 


g 


equation (B7) becomes: 


(26) 


^ <? | 

/P L +Py<P\ 

, COS 3 0 * 

\ - (2 - <p\' 

3(P ? - P v ) 3 1 

Vl - P v ) 

3 

> 

Q. 

1 

a 
J 


cos 6 * sin2 9 * cp / Pl ~ Pv A - Z ° Sm " P v <P 

2 \P Z - P v / 2 a/Bo 


P Z -P v 


sin 2 g* 
<^Bo 


V / Py 

/ '/§o \p ; - p v/ 


w 
__J 

7TW*f(0*) 


1 - 

2 cp 


(Sin 2 0* (B8) 
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Regrouping, the density ratios of equation (B8) 


2p 3^ <P ! p l * p v <P \, 1 1 _ 2 P a <P - P t f - P y 9 2 + p t - P v 1 

3(p t - P v ) 3 \P l -P v ) 3 3 3 <P Z - P v ) 3 


[’ 


2<pp s - (1 +(p )pJ 


+ P t ( 1 - <p) 


3(p, 


P v ) 


- (B9) 

3 


And finally substituting equation (B9) into (B8), ? equation (B8) reduces to 


2W S - (1 + <P 2 )P V 

+ Pj(l - <p) ^ i cos S | n 2 

P t - (2 - cp 2 )p y 

3(P Z - P v ) 3 2 

_ Pi ~ p v - 


cos^ 0* 


<p 






+ 


V / p v \ 



7TW*f(0*) 



0 * 


(BIO) 
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APPENDIX C 

SURFACE TENSION SUPPORT (ISOTHERMAL-SPHERE) 

If a load is placed on an axisymmetric shell or membrane structure, 
it is in effect transmitted to a contact circle., In our case this contact cir- 
cle is represented by the liquid- vapor-solid, or fluid A-fluid B-solid inter- 
face. The surface forces acting at this interface are illustrated in figure 3. 

While surface tension, a y , and the wetting angle a characterize the 
interface, it is useful to think of surface forces in vector terminology. 

The components of surface tension acting to support the sphere and 
those tending to pull the sphere into the liquid are Illustrated in figure 5. If 
we define the adhesive tension at the contact circle to be which acts 

normal to the surface, then a balance of forces requires that 


a 


adh 


sin a 


(Cl) 


As the horizontal forces are balanced because of symmetry, vertical up- 
ward force, F , balances the vertical downward forces: 

0 = F a - 2 ttR 0 sin 0*(-a a ^ cos 0* + ct sin 0* - cr gv sin 9 *) (;C2) 

where 2ttR sin 9* is the circumference of the contact circle. But the 
o 

Young-Dupre equation gives 

a i V C0S “ + a sl = °sv (C3) 

and upon substitution of equations (Cl) and (C3), equation (C2) becomes 

F^ = -2tiR 0 ct^ v sin 0* sin (0* + a) (C4) 


Equation (C4) also results from assuming that cr^ y and a completely 
characterize the interface. Resolving a , the vertical component becomes, 


v cos (9* + a) - 


3 n 


- - a, y sin (0* + a) 


(C5) 


and integrating equation (C5) over the contact circle gives equation (C4). 
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APPENDIX D 

ESTIMATING 0* FOR NONOPTIMUM FLOATER 

Enter figure 4(a) with the density ratio (p e - p )/(p - p ) and or, 

L s v t vj | 


determine the critical Bond number Bo 


cr, 


If we assume that the Bond 


number change is effected by a change in R , then estimate a new density 
ratio 



(Dl) 


(P l ~ P v )g| 


Re-enter figure 4 with (p - p )/(p - p ) and a and determine a new 

L s v i v jo 


critical Bond number Bo 


cr. 


2 

Using Bo find 6*, and z n from fig- 

V I 


o 


ures 4(b) and (c), respectively. The technique appears to be valid for 
small changes in Bond number. 
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TABLE I. - SOME EXPECTED ASYMPTOTIC 
BEHAVIOR OF THE CRITICAL 
LEVITATION LOCH 


Bond number 

Density ratio and 0* 

Bo “*■ e ^ 0 

a — u 

P s ■ p v _ 3 
- Py 2Bo 

0 * - £ 

2 

2 

1 

p s " p v _ 3 

P Z ' p v 4Bo 

0* _ 37T 
4 

Bo — °° 

p s ~ P v _ i 

p l ~ P V 
0* — 7T 

For all Bo 
a-* 0 

p s ' p v _ j 

p l ~ p v 

0* — 7T 
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w s + / P v cos w dAj -/ p 0 cos 0 dA 2 W $ + / P v cos w dAj - / P t cos 0 dA 2 
-/T re sln0dA 2 -O - f o L y sin(0 !;: + a)dl = 0 



(A) FORCES ACTING ON THE NON- (B) FORCES ACTING ON THE ISO- 
ISOTHERMAL FLOATING-SPHERE IN THERMAL FLOATING -SPHERE. 
FILM BOILING. 


Figure 1. - Free body force balance for levitated spheres. 



(A) THE NONISOTHERMAL FLOATING- (B) THE ISOTHERMAL FLOATING-SPHERE. 


SPHERE IN FILM BOILING. 


Figure 2. - Detailed force analysis of floating spheres. 





FLUID SPHERE 

a 

o FLOAT 
o' 

WATER 


a FLOAT 
O' NO FLOAT 

ETHYLENE GLYCOL 

7 r 
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7T 

^ FLOAT 
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7 r 
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CHLORIDE 


a FLOAT 

TEFLON ON WATER 

3n/5 

o NO FLOAT 

ISOTHERMAL 




BOND NUMBER B 0 


<PrPy>9 R o 

og c 


Figure 4(a). - A comparison of the critical levitation locus 
and various floating spheres on liquid nitrogen. (A 
Teflon sphere floating on water is included for comparison.) 


UJ 



Figure 4(b). - Submersion angle 0,^ for the optimum 
density ratio. 





pt SUBMERSION DEPTH FORM 
MAXIMUM DENSITY RATIO 



Figure 4(c). - Submersion depth Zq for the optimum density 
ratio. 



°ADH MUST BALANCE 
o LV sin(7r-a) = o LV sin a 


Figure 5. - Surface forces acting at the contact circle. 
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